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Abstract. We introduce a general class of statistical-mechanics models, taking 
values in an abelian group, which includes examples of both spin and gauge 
models, both ordered and disordered. 

The model is described by a set of "variables" and a set of "interactions". 
Each interaction is associated with a linear combination of variables; these are 
summarized in a matrix J. A Gibbs factor is associated to each variable (one- 
body term) and to each interaction. 

Then we introduce a duality transformation for systems in this class. The 
duality exchanges the abelian group with its dual, the Gibbs factors with their 
Fourier transforms, and the interactions with the variables. High (low) couplings 
in the interaction terms are mapped into low (high) couplings in the one-body 
terms. If the matrix J is interpreted as a vector representation of a matroid, 
duality exchanges the matroid with its dual. 

We discuss some physical examples. The main idea is to generalize known 
models up to eventually include randomness into the pattern of interaction. We 
introduce and study a random Gaussian Model, a random Potts-like model, and 
a random variant of discrete scalar QED. 

Although the classical procedure 'a la Kramers and Wannier does not extend 
in a natural way to such a wider class of systems, our weaker procedure applies 
to these models, too. We shortly describe the consequence of duality for each 
example. 

PACS:05.50.+q, 6160. Cn, 75.10.-b 
keywords: Duality, Matroid Theory, Potts model. 

1. Introduction of the problem 

In statistical mechanics a notion of duality appeared first as a relation between the 
high- and low-temperature expansions in the two-dimensional Ising model pp. In 
both expansions the partition function can be expressed as a sum on a gas of closed 
polymers: in the low-temperature picture, these polymers are the boundaries of the 
ferromagnetic domains, so we have an expansion in the deviations from the ordered 
regime; on the contrary, in the high-temperature picture the polymers are the terms 
of a Fourier expansion which expresses the deviation from a disordered regime of total 
decorrelation. So, in a certain sense, the duality relates energetic degrees of freedom to 
entropic degrees of freedom. Moreover, under the assumption of the existence of only 
one critical value for the temperature, it was possible to locate the phase-transition 
point before the exact solution of the model by Onsager. 
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More formally, a duality transformation can be defined directly on the partition 
function of the model by means of two independent duality transformations, that is: 

• The variables are located on the /c-cells of a 13 dimensional cell-complex (e.g. the 
0-cells, or sites, for a spin model, the 1-cells, or bonds, for a gauge model), more 
precisely the physical state is a fc-cochain taking values in an abelian group. 
The interactions are located on the (fc + l)-cells (respectively the bonds and the 
plaquettes in the previous examples). The "algebraic" duality exchanges variables 
and interactions, and the dual physical state is a (A; + l)-chain on the dual of the 
original group. 

• The new variables have a gauge redundancy, and must obey a set of constraints: 
in other words, they must be a closed (fc-l-l)-chain. The solution of the constraints 
leads to the "integration" of the (/c+l)-chain to a (/c+2)-chain. After this, the cell- 
complex on which the model is defined is transformed into its dual counterpart. 
This "geometrical" duality of cellular complexes pictorially exchanges /c-cells with 
{D — k)-ce\\s, and the physical state is now a, [D ~ k — 2)-cochain. 

This is the general strategy, explained in detail, for example, in the review part of |2], 
or in 

In this paper we wish to clarify the first ("algebraic") part of the duality 
transformation, by putting it in a more general context. We consider a rather 
general class of statistical-mechanical models in which the variables take values in an 
abelian group; each interaction term is associated to a particular linear combination 
of variables, but need not arise from any particular geometrical entity (such as a bond 
or a plaquette). By relaxing the requirement for a cell-complex structure, we can deal 
with a wider class of statistical models. We think that this more general notion of 
duality can be fruitful in the study of disordered systems. 

This notion of duality is motivated, in fact, by the clauses-variables duality that 
we observed in a previous investigation in the zero-temperature behaviour of a model 
of random satisfiability 0j, which has the structure of a disordered diluted interaction. 

The paper is organized as follows: 

• In section 2 we introduce some mathematical formalism, in particular about the 
notion of duality and Fourier transform for generic abelian groups. 

• In the sections from 3 to 8 we discuss our duality under general aspects: in section 
3 we describe the most general setting for the duality; in section 4 we specialize 
the model to a formulation which is manifestly suitable for statistical mechanics 
models, and in section 5 we rederive the duality in this frame; in section 6 we 
show how to reconduct an extension of the original definition to the primitive 
formulation of the problem, which is reminiscent of a set of gauge fields coupled to 
the interaction terms; in section 7 we describe a particular extra-feature appearing 
in the model when delta constraints appears; finally in section 8 we introduce a 
sufficient criterium for self-duality. 

• In sections from 9 to 13 we show the application of these general ideas to some 
concrete physical examples: first, in sections 9 and 10, a sort of random gaussian 
model, and a variant with random constraints; further, in sections 11 and 12, 
a generalization of Potts model extended to random structures of interactions, 
and a discussion on how to recover, inside our framework, the traditional results 
for regular lattices; finally, in section 13, we discuss a concrete example of the 
general setting of section 6, which is reminiscent of a scalar field theory coupled to 
electro-magnetism, both in the case of regular lattices and of disordered systems. 
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• In section 14 we outline some conclusions, and introduce the perspective for a 
further generalization of the procedure, using different mathematical structures 
instead of an abelian group structure and the Fourier transform. 

2. Dual for an abelian group, Fourier transform and group 
homomorphisms 

In this section we introduce the Fourier transform for locally compact abelian groups, 
by means of the theory of characters, following the textbook of Kirillov and Gvishiani 

First, we introduce the concept of dual group. Let G be a locally compact abelian 
group (we always write our abelian groups additively) . Let TZq be the set of irreducible 
representations of G. It is well known that the irreducible representations of G are 
all one-dimensional: thus, for each p G TZg, the representation matrix Up{g), is one- 
dimensional and coincides with the character Xpio)- The set TZq is itself an abelian 
group under ordinary composition (multiplication) ; we write this group additively and 
call it the dual group G. It is understood in the following that latin letters g, h. . . are 
elements of the original group G, and greek letters like p, r. . . are elements of the 
dual group G. If we adopt the notation 

Xp(g) - U,{g) = e'^-^"^^) , (1) 

the following bilinear relations hold: 

i9g {p + T, g) = i?G (/o, a) + (r, g) ; (2a) 
7?g(p, 9 + h) = i^aip, g) + ^cip, h) . (2b) 

It can be shown that the dual group of G, denoted G, is naturally isomorphic to G. 

In a consistent choice of notations, the characters of the dual group are derivable 
from the antisymmetric relation 

dG{p,g)^-^a(9,p), (3) 

which is equivalent to 

XsCp) = Xp(ff)- (4) 

Given a function / from G to the complex plane, its Fourier transform / is a function 
from the dual group G to the complex plane, defined by 



/: G^C /(g) = / e"'e(s.p)/(p) (5a) 

f:G^C /(P)= /e'''^('''^)/(5). (5b) 



where and J* stand for a properly normalized Haar measure on G and G (see 
below). Subscripts on the bilinear forms z^g and i^g will be omitted when clear. 

We introduce a notion of invariant (5-function compatible with the measure 
introduced above: 



/ Xpig') Xpia) = Saig' ~ g) ; Xp'ig) Xp{g) ^ 5^[p' - p) 

I p Jg 
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Table 1. Convention for the normalization of the Haar measure adopted in this 
paper for various pairs (G, G) of abelian groups. It is understood that for the 
direct product of groups the Haar measure is simply the product measure of the 
original ones. 



which have the natural property 

/ fig) Scig' g) - f{g') ; /* /(p) <5g(p' - p) = Kp') ■ 

Jg J p 

Inside the wide class of "physically natural" groups with the properties of being locally 
compact and compactly generated, the most general abelian group is the product of 
R, U{1), Z and Z^ groups, so the suggested choice of normalizations of TableQlcovers 
all the possible situations. | 

A homomorphism ip from an abelian group Gi to another abelian group G2 is 
a map which preserves the group operation ip{xi +Gi 2^2) = 'fiixi) +G2 vi^2), and 
maps the identity to the identity, Lp{0) — 0. We denote by Hom{Gi, G2) the set of all 
homomorphisms from Gi to G2. Note that a homomorphism must preserve additive 
inverses, i.e. g) = —ip{g). This follows from ip{g) + ip{—g) = <p(0) = 0. 

As we have seen in the previous paragraphs, from theory of characters we can 
introduce the dual groups Gi and G2. Following Kirillov ( 5 , Part III, statement 631), 
it is indeed possible to introduce a notion of dual homomorphism if'^ E Hom{G2, Gi) 

I This fact derives essentially from three general results 6 : 

• Every locally compact, compactly generated abelian group G is topologically isomorphic with 

X X K for some nonnegative integers a and b and some compact abelian group K. 

• Let K he a compact abelian group and let U he s. neighborhood of the identity in K. There is a 
closed subgroup H of G such that H <ZU and K/H is topologically isomorphic with U{1)'^ X F, 
where c is a nonnegative integer and _F is a finite abelian group. 

• Let be a finite abelian group. Then it is isomorphic to the group direct product of cyclic 
groups of prime power order [Kronecker decomposition theorem). 
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^ 

The dual homomorphism ip^ acts according to the formula 

Xp(¥'(5)) = Xff(¥'^(p)) = X>^r(p)(5) ; g^Gi; p e G2 , (6) 
from which we derive 

dGA^{x),y)=dGA^^iy),x)- (7) 

It is easy to verify that (p^ is a homomorphism from G2 to Gi, as a consequence of 
the fact that d{p,g) is a bilinear form §. 

3. The general formulation of the duality 

In this section we wish to introduce a notion of duality for a wide class of systems in 
classical statistical mechanics. 

In the most general form, given two abelian groups Gi and G2, we will consider 
the system defined by 

• a function Fi : Gi ^ C; 

• a function F2 : G2 ^ C; 

• a homomorphism G Hom{Gi, G2); 

and the Gibbs weight of a state a; G Gi is given by 



exp 



(8) 



The partition function is defined as 

Z 



FuFr,v- I Fi{x)F2{p{x)) . (9) 

Applying Fourier transform to the function F2 we obtain an algebraic relation: 

FAx)F2{^{x)) ^ I I FAx)F2{y)e'''^-^^^^^-^y^ 

= f f F2(y)i^i(a;)e'''«i(^"(^)'-)= / F2{y)F,{p^ {y)) , 

J y J X J y 

from which we deduce that 



(10) 



Note that in the dual model the roles of the group Gi and G2 are interchanged, and 
the homomorphism tp g Hom{Gi, G2) is replaced by p'^ E Hom{G2, Gi). 

§ The homomorphism property descends from the biUnearity properties 
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4. Definition of a model with many degrees of freedom 

In this section we specialize the procedure described above to a more physical frame, 
in which extensive parameters appear both in the choice of the groups Gi and G2 and 
of the functions Fi and F2. 

We introduce two extensive parameters: the first one, r, is called the number of 
variables, while the second one, (n — r), is the number of interactions (they were called 
clauses in 4 ). 

Given a group G, we adopt Gi = G®"" and G2 = G®^"-''). A state in Gi is 
specified by a vector x = (xi, . . . ,Xr), and the homomorphism ip is specified by a 
r X (n — r) matrix J, with elements Jik G Hom{G, G) acting on the left, such that 
{xJ)k means X]i=i Jik(xi). By convention we let the indices k run from r + 1 to n. 

The functions Fi and F2 will be chosen to be factorized on the entries of the 
vectors: 

r n 

Fi[x]^Xlf^{x,)- F2[y]= n h{Vk)- (11) 

i=l fe=r+l 

The weight functions /i , . . . fr are thus "one-body terms" (generalized magnetic 
fields), while the fr+i, ■ ■ ■ fn sue "interaction terms". 

Note that we have chosen the groups Gi and G2 to be direct products of the 
same group G. To avoid confusions in the application of the conventions of Tabled 
we specify that 



x]£Gi 



■dGi{x',x) ^2^dG{x'^,Xi) ; 

i 

f daAy' ,y) = -Y.'^G{yk,y'k) . 

■^[j/]eG2 •'Vr+i Jyn 1^ 

We will sometime use a notation in which the two functions Fi [x] and F2 [xJ] are 
"merged" into only one function F[xB] — Hj^i fj{i^B)j)- 

Denote with 1 the identity homomorphism ipi G Hom{G, G), such that ^pi{x) — x 
for any x, and with the null homomorphism ipo G Hom{G, G), such that 1^90(2;) = 
for any x. The r x n matrix B is thus defined as 

1 i = j; 1 < j < r; 
B^J = '{ ^ ^ j; 1 < J < r; (12) 
Jij r + 1 < j < n; 



which can be written simply as 



B 



Ir 



J 



(13) 



The invariance of the partition function under change of variables in the integration 
corresponds in terms of the matrix B to the application of a special orthogonal 
transformation from the left (change of basis) 

Z[B] = [ F[xB] ; Z[B] = Z[OB] V O G SO{r) . (14) 
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It results that, beyond the particular choice of the matrix B, a more fundamental 
notion of "pattern of interaction" can be established, which is independent from the 
choice of basis. This is the analogous of what happens in linear algebra, where the 
notion of linear independence is intrinsic with respect to the choice of basis. The 
introduction of this concept in mathematics dates up to the paper of Whitney |7], 
and is the foundation of a theory called Matroid Theory (cfr. [H| for a textbook). 
In particular, matroid is the name given to the intrinsic mathematical structure 
corresponding to our pattern of interaction, and a choice for the matrix B is called a 
vectorial representation for the matroid. In the case in which the first r x r block of 
the matrix is the identity, we say that _B is a standard vectorial representation for the 
matroid. Note that our particular choice puts B automatically in standard vectorial 
form. 

In this frame, we can develop a physical intuition on the structure of our model. 
We could consider the r entries of the vector x as the variables of a statistical mechanic 
system. The first r functions fi{xi) collect the one-body contributions to the Gibbs 
factor, like one-body measure terms, or magnetic fields, or chemical potentials. The 
last {n — r) functions fk{{xJ)k) acts like {n — r) terms of interaction, encoding 
informations like the strength of the couplings, or the nature of the clauses, etc. 

The matrix J describes how the variables interact. In the more general case, 
the entries of J are homomorphisms from G to itself. A specialization of this case 
is the one in which the entries of J are in a commutative ring with unit, TZ, and 
G is an 7?,-module. A further specialization is the case in which G is an arbitrary 
abelian group and 7?. = Z. In fact, every additive group is an 7?.-modulc with Z, where 
'n • g = g -\- ■ ■ ■ + g (n times) . 

An example of this case is a scalar field theory in which the r variables (pi take a 
real value for each site i of the lattice, and the interaction terms depend on [pi — 0j)^, 
where (i, j) is a bond of the lattice. Another example is a Gauge Potts Model, in 
which the r variables Si take a value in Zg for each bond i of the lattice, and the 
interaction terms depend on the linear combinations X^ieOp where p is a plaquette 
of the lattice. In fact in both these models all the coefficients of J are in {0, ±1}. In 
the first case, the entry Jy is non-zero when the site labeled by i is an extremum of 
the edge labeled by j, while in the second case this happens when the edge labeled 
by i is on the border of the plaquette labeled by j , and the sign is determined by the 
choice of orientation for edges and plaquettes. 

Nonetheless, other mathematical structures can be considered. A case which is 
not included in the situation above (G is an 7?.-module) is when the entries of J are 
in a field F, and G is a finite-dimensional vector space over F (say, F'). The field F 
could be, for example, the real field R, the complex field C, or a finite field GF{q^). 

Moreover, we stress the fact that the generality of our formulation does not require 
any underlying geometrical structure for the system, and applies also to situations, 
typical of disordered mean-field models, in which such a structure does not exist. 

5. Duality for product groups 

Let us now repeat the procedure of the previous section specialized to our case, and 
write the corresponding duality transformation. The Boltzmann weight is 



n 

exp(-7^G;/;s(.x)) = n/^((^^)j) 



(15) 
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The partition function is 



Zgj^b = / exp{-nG:f;B{x)) ; 

J[x] 

and, applying Fourier transform, can be restated as 

J\x] , ,1 J\x\J\y] 



(16) 



fe=r+l 



k=r+l 



n r n r 

n A n (^') = / n (y^^ n ^ ((-2^"^"^) 



'fel-'Fl k=r+l 

from which we deduce that 

where the dual i? of the matrix B is intended exactly as in matroid theory (see 
chapt. 2). 



B 



J 



B 



(18) 



In the most general case of Jy S Hom{G,G), the quantities {J'^)ji are the dual 
homomorphisms Jy defined according to © . 

The duality interchanges the role of variables and interactions. Furthermore, if 
the functions fj are almost constant (low coupling, or high temperature), the functions 
fj are very peaked (high coupling, or low temperature), so the duality relates a regime 
of high (resp. low) temperature in the one-body terms with a regime of low (resp. high) 
temperature in the interaction terms, and vice versa. 

6. The coupling with gauge fields 

In this section we slightly generalize the formulation of the problem introduced in 
sections 13 and ^ and then we show that, with a change of variables, it can be recast 
in the original framework. 

This analysis has one main reason: in section El we will describe the physical 
examples both of a scalar field theory, and of a discretized free electrodynamic theory, 
and show that our duality holds for both systems. Then we show that the duality 
holds also for the system in which electrodynamic is coupled to the scalar field in the 
traditional gauge-covariant way. 

We want to show that this is a particular case of a general recipe, introduced in 
this section, and valid for the cases in which we introduce a set of gauge fields coupled 
to the terms of interaction. The change of variables corresponds to convert the basis 
of scalar fields and gauge fields into the basis of fields and covariant derivatives. 

Given 2k abehan groups Gi, . . . , Gk^ G'^, . . . , G'j,, we will consider the system 
defined by 

• k functions fi from the group Gi to the field C; 

• k functions 'cji from the group G[ to the field C; 

• 2fc + 1 homomorphisms ipi between the groups Gi and G'{. 



0- 



V-k 



G, 



V-(k-l) 



■G[ 



Gi 







(19) 
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and the Gibbs weight of a state x — {xi, . . . , Xk) E is 

fe fe 
exp [-H/,g;y,(a;)] Y[ ^ f-i{xi+i))Y[9j{fj-i ■ ■ ■ (fioixi)) ; (20) 
i=i j=i 

(it is understood that Xk+i — and (p-.k{xk+i) = 0). 

Consider the change of variables Vi = Xi — (/9_i(a;j+i), with Jacobian 1. If we 
introduce the k x k matrix of homomorphisms 

Ipij = (p_(i_i)(^_(i_2) • • • (fij-^l , (21) 

where the homomorphism V'ij goes from Gi to G^-, we could restate the Gibbs weight 
above as 

exp[-nfrg;^M] = (2^^)11 5. (E^y(y')) ■ (22) 

As a matrix of homomorphisms from the set of groups {Gi} to the set of groups 
{G'j} is a homomorphism between the product groups Gi and G'j, we are in the 
conditions to apply of the general procedure we have outlined in section |3| 

The dual model interchanges the groups Gi with the groups G^, the functions fi 
with the functions gi. As the transposed homomorphisms are given by 

i^^^h = ii^J^f = ■ • • : (23) 

it corresponds to interchange the homomorphism ipi with the transposal homomor- 
phism iy9^j. In this specialization of the general procedure, the relation 1)1 0|l would 
read: 

Zf^Mv.} = ^gj-W^j- (24) 

In the case we deal with product groups, as described in section^ the matrix B will 
not be still in the standard form representation of Matroid Theory when expressed in 
the X basis. Nevertheless, it will be in a clear block form with k rows and 2k columns 
of blocks. If we call Ji the linear applications corresponding to the homomorphism 
ipi, the block Bij would be given by 



B. = { . :zir'r:i_ (25) 



The matrix B can be easily diagonalized. This results in the change of basis from 
X variables to y variables in the sense described above. The resulting matrix B' will 
have blocks B-^ given by 

{/ i = j < k ' 

^^J^l<k; (26) 

J-{i-i)J~{i-2) ■ ■ ■ Jj-k-1 j > k . 

Now that B' is in the standard form representation, we can apply duality, and then, 
inverting the diagonalization procedure described above, restate the dual problem in 
the original block form. The transformation interchanges the matrices Ji with J^j, 
with a minus sign on Jq, which goes into —Jq- 
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7. The solution of delta constraints 

A general feature of our duality in the factorized case is that the resolution of 
constraints taking the form of delta functions has a dual counterpart in the removal of 
trivial constant factors, and vice versa. This feature can be seen in a transparent way 
in a simple unstructured problem, shown in section [TUI but it is more commonly used 
when the duality is combined with other geometric properties of the systems, as, for 
example, in the case of models defined over regular finite-dimensional cell-complexes. 
Such a situation will be described in section IT^ 

Consider a system with r variables and n — r interactions, in the standard form 
described above, in the case in which the first r functions fi{x) are constant. For 
simplicity, say that the last (n — r) functions fkiv) ai'e all equal to a certain f(jj). The 
model is described by the Gibbs weight 

n 

eM-n{x)]= n 9{{xJ),). (27) 

j=r+l 

Following the duality procedure, the dual Gibbs weight is given by 

r n 

e^p[-H{0]=l[S{{-^J^),) n m- (28) 

i—l j—r+1 

In the following, assume for simplicity that the incidence matrix J has maximal rank. 
X In this case it can be written in standard form J = (/, J') by mean of a change of 
variables, and a relabeling of the last (n — r) indices. So we can rewrite the Gibbs 
weights above as 

2r n 

eM-'H{x)]= n /(^.) n /((^'^').); (29) 

j—r-^1 J— 2r+l 

r 2r n 

eM~'Hm = X{5{-^T''-{^^J'^)^) n /(^f') 11 /(€"'); (3o) 

i=l j=r+l j=2r+l 

where in the last expression we already stressed that the variables ^'"'^ form an 
independent set, while the variables ^'^"f can be expressed in terms of them when 
solving the delta constraints. In this case we obtain 

2r n 

exphw(o]= n fii-r^j'^w n /(^i"')- (31) 

j=r+l j=2r+l 

This is exactly the Gibbs weight obtained applying duality directly to the weight in 
the (|29(l . considered as a system with r variables, labeled with the indices from ?- + 1 
to 2r, and n ~ 2r interactions, labeled with the indices from 2r + 1 to n. 

I Matrices of non-maximal rank r' < r would only imply trivial overall factors of the kind |F| \ 
or the necessity of fixing "by hands" (r — r') variables, as for example in planar Potts model, where 
a factor q originates from this fact, or in the XY-Model - Coulomb Gas relation, in which we should 
avoid to sum over the variable associated to one special site, in order to fix a "ground" for the 
potential. 
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t 

(Duality ) 
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' 
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I, J 
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f 

(Duality) 
} 

^JT. I 
f 



8. A sufficient condition for self-duality 

Duality is particularly useful whenever it is possible to recover that a system, under 
suitable conditions, is self-dual. 

In the following, let us denote by X the instance of a problem (say, the lattice, 
the set of couplings,. . . ), and by cr a configuration of the system (say, the values of all 
the spins in the system). 

Duality means that there is an involution {X} {X} such that 

Z[X] = f[X]Z[X] , (32) 

with f[X] an analytic function of X. 

Of course a sufficient condition for an instance X* to be self-dual is that it is a 
fixed point of the involution 

X^ = X^ . (33) 
If there is a group of invariances U such that 

VueU Z[uX] = Z[X] , (34) 

then a weaker condition for self-duality can be stated. It is sufficient that some u &U 
exists such that 

uX^ = X* . (35) 

At this point, we make a short digression on disordered systems. For these systems the 
variable X describing the instance is itself a random variable, distributed over some 
space {X} with some normalized measure fi{X) given a priori. We are interested in 
the thermodynamic quantities (that is, thermalized over a) of a system in a given 
instance X, furtherly averaged over the distribution of disorder /Lt(X). For example, 
the average free energy, up to /3 factors, is given by 



F = logZ = j diJL{X) log Z[X] . 



Source variables can be introduced to calculate the other thermodynamic quantities. 
For example, to calculate the average magnetization we can introduce a magnetic field 
h, and obtain 
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similarly, to calculate the average two-point function {(Jidj) we can introduce a vector 
of local magnetic fields {hi}, one per site, and obtain 

dhidhj 



{a.,a,) - / MX) = ^^ri^F{h) 



h=0 

The action of duality on the space of instances {X} naturally induces an action on 
the space of normalized measures over {X} 

^ ^ ^ ^ OX , , ^ , , 

oX 

The condition for self-duality now applies to measures. The trivial condition for a 
measure IJ*{X) to define a self-dual problem is 

9{x) = ^l*{x). 

Analogously to what we did above, if we have a group hi of invariances for the problem, 
we can make this condition weaker. If we have a normalized Haar measure on the 
group, the condition which exploits the symmetries of the problem is 



^i*{uX) = / ^i*{uX). 
u Ju 

Of course this condition can be fulfilled, by using the invariance of the Haar measure, 
if for any fixed X there is an element of the group U such that 

1?{uX)^ ^i*{uvxX). (36) 

In practice, in the case of a disordered system, we have a measure on the space of 
n-uples of functions {/i}, and of r x (n — r) incidence matrices J. In the case in 
which r = n/2, all the first r functions fi are equal to a certain f{x) and all the 
last r functions are equal to its Fourier transform f{x), the dual set of function is 
identical to the original one, and we only need to check self-duality on the measure 
for the incidence matrix, /i(J). 

Let the functions f{x) and /(x) be invariant under a subgroup of the 
automorphisms of the space G, ^ C Aut(G'). 

The elements u of the group of invariance lA are such that Zf -j — Zf-^j- 
They are of the form u — (tt; {oi}, {0^}) with tt e S{i,....r} x &{r+i,...,n}, where 
©A is the permutation group on the set A, and a^, Ok in A, with i — 1, . . . , r and 
k = r -\- 1, . . . ,n. The action of u on a matrix J is 

from which we have that a consistent group composition is 

u'u = (tt'tt; {a^/(j)a-}, {aj,a^/(fc)}) . 

The condition for self-duality (|36|l above, in our case can be restated as the existence 
for each matrix J of an element uj U such that 

^i*{uJ{~J^j)^^,*{J). (37) 

Note that the Jacobian of the transform J — J-^ is one. Exploiting the group U 
of invariances for the problem, we can restrict our attention to measures which are 
central functions for the group U. If we introduce the two-element group generated 
by the transformation l : J ^ —J^, the self-dual measures are the (fewer) central 
functions on the (larger) group U x: {e, t}. 
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9. An easy example: the Gaussian model 

As a first easy application of formula H17|l we will consider a simple model exactly 
solvable via Gaussian integration. Consider the system described by equation H15|l . 
specialized to the case in which G = M and the functions {fj} are 

=e-^'/2<T= i = l,...,r; (38a) 

/,+fe(a;) =e-=^'/2r2 fc = l,...,n-r. (38b) 

We will introduce a notation ad hoc to label the partition function of this system: 

ZGa.uss{J;{o-i}, {Tfc}) = ZR B,{fj} , (39) 

where the expressions (|38|) above for the {/,} and the vectorial representation 
B = [Ir\J] are understood. 

We could rewrite equation p5(l in the more familiar form 

n 

eM-nc-jM^)) - n m^B),) = e-i^'{s-^+JT-^Jn^ , (40) 

where Sw = (JiSui and T^fe' — TkSkk'- Using the convention of Tablenfor the measure 
in the case A = M (that is, writing down the proper number of (27r) factors), and 
applying the formula for Gaussian integration, we have for the partition function 

J i/det {Ir + M^M) 

where we defined 

M = T-^J^S. (42) 
From equation (|17(l we have that the dual partition function is 

^Gauss(-J^;K7l},K-l}) - ^^ToT4ir' ^Gauss(J;{aJ,K}), (43) 

(27r) 2 det 5 

while, applying the formula for Gaussian integration on the dual formulation, 

(44) 

So the consequences of the duality can be restated in the algebraic equation 

det {Ir + M^M) = det {in-r + MM^) , (45) 
which is indeed true, although not totally trivial. § 

§ Consider the formal power series in z: 



ln(det(/ - zM^M)) = tr{ln{/ - zM'^M)) = - V — tr(Af'^Af )" 

n 

n— 1 

= - — tr(MM^)" = tr{ln(/ - zMM'^)) = ln(det(/ - zMM'^)) . 
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10. A randomly constrained Gaussian Model 

In this section we introduce a model in which some Gaussian distributed variables in 
are subjected to a set of random constraints. 

This model is a first example in which the matrix of the homomorphisms B really 
contains non-trivial group homomorphisms. Furthermore, it exhibits in a clear way the 
general feature of our duality described in section|7| the resolution of delta constraints 
has a dual counterpart in the removal of trivial constant factors. 

The model is defined as follows. Let Xi, i — l,...,r, be ZJ-dimensional real 
vectors, and introduce the (n — r) linear combinations 

2/." = EE«" (46) 

i=i 0=1 

with = r + 1, . . . , n, defined by the matrices Jik- We will consider the case 2r > n. 
The Gibbs weight of a configuration is 

da;exp[-Hj(a;)]-d^a;i---d^a;, exp(-— ^x^) • J| %fc,0). (47) 

1=1 k=r+l 

Working in the spirit of disordered systems, we would have a probability measure 
over the quenched disorder variables Jik ■ Introducing the Haar measure d/i d over the 
group of special orthogonal transformations SO{D), a possible choice is to take the 
variables Jik as i.i.d. variables, distributed with d/i/j: 

df^iJ)= n ^^^D{J^k). (48) 

'i=l,...,r 
fc=r-t-l,...,n 

The problem has a large gauge group, which allows for the reduction of the matrix 
J in a standard vectorial form. We have mainly three class of manipulations on B 
which do not modify the value of the partition function: 

• permutations of the columns: 

Jik ^ Jia{k) ^ ^ r 5 

• global rotations on rows and columns: 

4" KU'^'^'^i^)' ^U) e SO{D), l<j<n; 

• full Gauss algorithm on the rows: 

J^k ^ QjJjk CeSOir). 
Thus, by change of variables we can put the matrix J in the form 




Remark that this can be done always when the rank of J is maximal, and this happens 
with probability 1 for our measure dfi{J). 
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This procedure of diagonalization corresponds to the algebraic resolution of the 
delta functions. In fact, the linear combinations are 

r D 2r-n D 

iy')k = E E ^'(J')'k = E E - ' (49) 

i=l 13=1 i=l (3=1 

and, as the matrix C above is orthogonal, the linear transformation which relates the 
new {y'l,} to the old {yk} has Jacobian equal to 1. Defining 

2r—n 

xt^in-r)ixi,---,X2r^n-Y)^ ^ xfyf^"" , k = r + l,...,n, (50) 

i=l 

and a transposal matrix such that (i^"^)^f = Y-^", the Gibbs weight (|47|l can be 
rewritten as 

^ 2r—n r 

exp[~Hj(a;)] =exp[-— ( ^ ^ |x,(x;F)p 

i—1 7 — 2r— n+1 

51 

2r-n ^ ' 



exp 
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which is of the form of the "simple" Gaussian model, equation l|4()(l . with 2r — n 
variables, n — r constraints, and all the Ui and the equal to tr, and the measure of 
integration dx stands for d^xi ■ ■ ■ d^X2r-n- So, with a slight abuse of notations, we 
can write the partition function of this model as 

ZG.nss{Y, a) = j dx exp ( - i^x^U + YY^)x) , (52) 

We could apply the duality to this formulation of the problem, using the general result 
of equation obtaining 

ZG.nss{Y, a) = ^^-^ ' ZGaus.(-r^, a'^) . (53) 

We could also apply the general duality directly to the original formulation. In 
direct formulation we have r variables xi, . . . ,Xr^ with one-body functions ,fi{xi) = 
e-^i/(2<^ and (n — r) interactions, encoded in the delta functions gkiuk) = ^iUk^O)- 
So that the matroid B has the form B = {Ir, J). 

Accordingly to the general theory, the dual problem has {n — r) variables 
r]i, . . . ,r]n-r, with constant one-body functions gkivk) — S{r]k) = l/(27r), and r 
interactions given by 



fM) = a/I exp (-^) ; ef = E E ^^(-J^)t^ 

^ k=r+la=l 



The dual matroid is given by _B = (— J"^,/„_r), and the definition of the transposal 
matrix is (— J"^)^f = —Jik- this point, we can remove the identity block from the 
matroid, as the functions associated to its columns are trivial constants. Then we can 
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diagonalize the matrix —J^ using the "transposal" of the procedure described above 
for J (that is, moves on rows and columns are replaced by the correspondent moves 
on columns and rows), and obtain again the partition function on the right side of 
(ESI- 

We remark that, in this case in which the entries of the matroid B are non-trivial 
homomorphisms, we can check directly that, accordingly with the general recipe of 
section the entries Bki of the dual matroid are related to the homomorphism- 
transposal of the entries Bik (that is, the internal indices a and f3 are exchanged). 
Moreover, in agreement with the general description of section [7| we see how the 
resolution of delta constraints and the removal of constant factors are dual procedures. 



11. A generalized Potts-like model 

In this section we study the effects of the duality on a particular model which slightly 
generalizes the one we introduced in 0]. We deal with r variables Xi e {0, . . . ,q ~ 1} 
and a general hamiltonian with (n — r) Potts-like interaction terms provided by a 
matrix J 

n — r r 

-?^(a.;J,{Xfc},{M) = ^Xfc,5(^a;,J,fc) +^M(x.), (54) 



k=l i 1=1 

where the sums and the delta functions are intended on Zg, that is a; = a; + nq. The 
weights of the terms in the Hamiltonian are more conveniently expressed in terms of 
a vector v G M", defined as 

Vi{x) = e'^^ — 1 for i — l,...,r; 

Vr+kix) — e^'' — 1 for k = l,...,n — r; 

and of the functions 

j,{x)^\ + v,E{x). (55) 

If we introduce the matroid B = (/^ | J) , we can reconduct the problem in the general 
form (fT5|l : 

n 

exp(-7^z„B,{/.}(^)) ^Ylil + vM^Bh)] . (56) 
We will introduce a notation ad hoc also for the partition function of this system: 

Zpous{J;{vj}) = Zj^^^B,{f,} , (57) 

where the expressions H55(l above for the {fj} are understood. 
From the form of the Fourier transform of f{x) = 1 + vS(x): 



/(O = - (1 + j . (58) 

and with the convention of Tabled for the measure in the case X = Zg, we have 
the relation between the partition functions 

^Potts(- J^; {q/v,}) = --ni^i^Zpott.(J; {vj}) . (59) 
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Let us specialize the general duality transformation we have obtained to some simple 
cases, related to the model that we introduced in Ref. j^. First of all, we will consider 
the case in which the magnetic field is constant {vj = w for all j = 1, . . . , r) and all 
the couplings are equals {vj = w for all j = r + 1, . . . , n). The equation H59|l in this 
case gives 

ZpoUs{-j'^; q/w, q/v) = ^^^^^„„^ ^Potts(J ■,v,w). (60) 

We can give a simple mathematical interpretation to the limit t; ^ 0, w — > oo of the 
formula above. In the two dual models, the partition functions, properly regularized, 
just count the number of solutions of the linear equations xJ = and y{—j"'") = 0. 
In this limit case the duality relation states that the rank of the matrix J is equal to 
the rank of — J-^, which is a trivial fact of elementary geometry. This model is exactly 
solvable for the ensemble of random matrices on a finite field, that is we know the 
exact probability distribution for the free energy and other physical quantities, and 
the duality relation plays indeed a role in the solution of the problem jl]. 

A particular case of (|60|) is when only v ~f oo (and then w' —^ 0) where again the 
partition functions should be properly regularized. The explicit expression is 

5:5(J.)(l+u;rS.(i-^(.0)^^]^ 

X 7) k 

therefore, the duality relates a sum, restricted to the set of solutions of a linear 
system Jx — 0, where the difficulty of the problem is contained on the external 
field contribution, to a sum, with no external field, on all the states of the spectrum, 
where the difficulty is contained in the non-trivial correlations between the degeneracy 
of the energy levels. A quantitative study of this system will be done in a forthcoming 
paper 

We can use the general strategy of section |H1 to investigate some sufficient 
conditions on the measure /xo(J) for equation to relate the partition functions of 
the same model at different values of the parameters v and w. 

As the weight functions involved are of the form "delta + constant", with no 
disorder parameters, the subgroup of automorphisms A C Aut(Zq) such that for a 
a G Awe have f{a{x)) — f{x) is isomorphic to the classical group Z*. || That is, if a 
is an integer coprimc with q, the corresponding automorphism on acts as a{x) — ax 
mod q. 

The group of invariances obtained is sufficiently large that, following the strategy 
described in section |51 the self-duality condition can be verified for almost every 
reasonable ensemble of matrices. For example, this is the case of random matrices, of 
sparse matrices, both with non-zero elements chosen randomly among {1, . . . , q— 1} or 
all equal to 1, and also of matrices with non-trivial correlations between the entries, 
if the correlations are symmetric under the exchanging of rows and columns, like 
the ensemble of incidence matrices of random bipartite graphs with vertices of fixed 
connectivity. An exception is the standard ensemble for the study of XOR-SAT 
problems ^OIj which is essentially a particular case of our model, but the matrices J 
are constrained to have a fixed number of non-zero entries per row, and do not have 
any constraint on the columns. 

II The group Z* is defined as the integers in {1, . . . ,q — 1} which are coprime with q; the group 
operation is the product, and the standard modulo q identification is understood. 



(61) 
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We note that, in the case of self-dual measure on the ensemble of matrices J, and 
V — w, the self-dual coupling occurs for v = ^/q, which is indeed the critical point 
also for the ordinary Potts model on a geometrically self-dual planar graph. We will 
show in the next section how to recover explicitly the self-duality in this case, in which 
other "geometrical" manipulations are required. 



12. Potts model on graphs and on planar graphs 

The model discussed in the previous section is a generalization of ordinary Potts 
model. In this section we will show how our duality relation can be reconducted to 
the common duality of Potts model on planar graphs 1 1 . 

Consider a connected graph A, at this stage not necessarily planar, with sites 
s g 5, bonds b £ B and elementary cycles (or loops) I G C. Say that the number 
of sites, bonds and loops are respectively S, B and L. If we consider the graph as 
embedded on an orientable surface, the set C of the loops could be chosen to be the set 
of the plaquettes ("local" loops), plus the set of the generators of the first homology 
group, that is, loops which turn around its handles. We remark that, at this point, 
this division is not necessary. 

We will give an arbitrary orientation to the bonds and to the loops. That is, for 
each bond b, Sin(6) and Sout(&) will be respectively the first and second extremum of 
the bond. Furthermore, each loop I will be associated with a i3-dimensional vector 
Jb{l), such that, having arbitrarily chosen an orientation for the loop, Jb{l) is the 
number of times that the loop I is occupied by the bond b in the same orientation, 
minus the times it is occupied by b in opposite orientation. 

The Potts Model on the graph A is defined by the set of variables ag £ 
{0, . . . , (7 — 1} on the sites, and a set of interactions with coupling constants Kj, on the 
bonds. % The Hamiltonian is 

-n Potts {cr; A, {Kb}) = Vk,,(5(ct^._^(;,) -cr5„„t(b)). (62) 

Graph f — f 
beA 

If we introduce the functions 

fsix)^l; (63a) 
fb{x) ^1 + VbSix); = e^^-l; (63b) 

and the matroid B — {Is\Jsb)^ with J a S x B matrix defined as 

{1 s = Sin (6) ; 
-1 s-Sout(&); (64) 
e.w. ; 

the partition function 

(T s b s 

^ We introduce an additive group structure Z,q for the q "colours" of the model, although ordinary 
Potts model Hamiltonian has a &q total permutational symmetry between the colours. We stress 
again the important fact that in general the symmetry group of the model is not related to the group 
introduced to perform the Fourier transform. 
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is in the general form H16|l . 

Now we remove the useless identity fiinctions fs{x), and make the change of 
variables 

Xb{a) = (Ts.„(fc) - crs„„,(h) . (66) 

The weight of each configuration a is now just Y\^^ fb{xb{<y))- Next, we want to perform 
the sum of the partition function on the new variables x, instead that on the old a. 
To this aim, we need to know how many configurations a correspond to a certain 
configuration x. The answer is simple: they are exactly q if the circuitation of x on 
each loop is zero (that is, Jb{l)xb = mod g), + and zero if the circuitation is 
non-zero for some loop. So we have 

o- X leC b 

and, with the definitions 

fbix)^l + VbS{x); fi{x)=d{x); (68) 

Jbl = Ml); B' =(Ib Jbl); (69) 



we have again a system in the general form H16|l . At this point we can perform the 
duality transformation. According to the general formula (|17|l . and to the expressions 
for the Fourier transforms of fb{x) and fi{x), respectively 

M) = -(^ + -m); M)^^; B'={-JibiL); (70) 



q \ Vb 

the result of the transformation is a system in the original general form 1()5|) . If the 
matrix — Jib is of the kind of Jgb, with exactly two non-zero entries per column, a -1-1 
and a —1, we can reconstruct a dual graphical version of the original Potts model. As 
shown in detail in the [Appendix A| this always happens when the original graph is 
planar. In a few words, when a graph is planar, every bond appears in exactly two 
loops (the plaquettes on the two sides). If we choose a uniform orientation for the loops 
(all clockwise or all counter-clockwise), it appears in the two loops with opposite sign. 
So, in this dual model, two neighbouring plaquettes of the original graph appear to 
interact with the dual coupling constant of the common bond («[, — q/vb). If we want 
to reconstruct a graph such that the system has the variables on the sites, and the 
interactions on the bonds, the desired graph is exactly the dual graph of the original 
one, in the sense of cell- complexes duality. 

The procedure of this section can be applied with minor modifications to obtain 
the Xy-model - Coulomb Gas relation {12i {131, provided that the XY-mode\ 
Hamiltonian is in the Villain form. Essentially, the only difference is that the couple 
of dual groups is (Zq,Zq) for the Potts Model, and (C/(1),Z) for the XY-mode\ - 
Coulomb Gas relation. 

13. The coupling with gauge fields 

In this section we describe the application of the procedure described in sectional to a 
very well known physical situation. The abstract framework that we introduced should 

+ More generally, for A not connected, we have a factor q for each connected component. 
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now assume a more transparent meaning: that procedure mimics what happens if we 
introduce a new set of gauge fields, coupled with the terms of interaction, and then a 
new term of the Hamiltonian depending from the physical content of these new fields. 

The paradigm of this procedure is the introduction of minimal electrodynamic 
interaction in a scalar continuous field theory, with lagrangian 

£0 = Jd<j>[V{cl))+d^cl>d^<l>] . (71) 

The derivative operator has to be modified into a covariant derivative operator 
Dfj,, defined by the action on the fields as -D^0 = dfi(j) — -A^. Then we have to 
add the dynamics of the gauge field, Ffj^^F^", where the field tensor is defined as 
Ffj.i' = df^Aiy — di,A^. The new lagrangian will be 



-Cqed = ld(j)dA 



(72) 



If we consider the statistical mechanic model whose Hamiltonian derives from the 
lattice discretized version of this model (as it is done, for example, in we replace 
partial derivatives 9^(/)(x) with finite difference operators {d4>)n-p, = 4'n+'ii ~ 4>m so 
that, as usual, the interactions lie on the bonds of the lattice. The linear combinations 
{d4>)n-fi are the arguments of the interaction terms in the Gibbs factor. 

The discretized gauge fields A„;^ couple to the terms of interactions, in such a 
way that the finite difference of the fields, (90)„;^, are modified into the covariant 
derivative (£)(/))„;p = {dcj))n-^i - An-^. 

The discretized field tensor Fn-^^i, = + An+ji-i, — Anj^^.^^ — An-^ takes the 
circuitation of the vector field An-^ around each elementary plaquette (n; iJ,v) of the 
lattice, and thus is a linear combination of the entries A„;^. 

The current quadrivector is given by the divergence of the electromagnetic tensor, 
jn,fj. = J2,yiFn,fj.L' — Fn-v,fiv), and its divergence, which is automatically zero, is given 
by (V • j)„ = Y.t,{in4^ - in-t,,t,)- 

We can generalize this frame to an abstraction of the geometric structure of a 
lattice, in which we have a set of sites {s}, a set of bonds {6}, a set of plaquettes {p} 
and a set of dual sites {s'}. In the previous case we would have for the sites {s} — {n}, 
for the bonds {b} = {(n, /i)}, for the plaquettes {p} = {(n, jJLv)}, and for the dual sites 
{s'} = {{n,^ivp)}. 

The action of gradient and divergence operators is now performed via some 
adjacence matrices. In particular we will introduce three matrices J^f,, Jbp and Jps', 
such that, given the two sets respectively of scalar fields 0s and of gauge fields A^, we 
have 

(D0)b = (j)sJsb ~Ai; Fp = AbJbp ; (V • j%' = FpJps' . (73) 

We will be interested in the theory with Gibbs weight 

exp{-Hc^M^,A)) =]JU4>s)l[MD<l>)b)l[gp{Fp)Y[gA{^ ■ j)s') , (74) 

s b p s' 

with the four set of functions given by 

fs{x) = 1 ; fb{x) = exp - ^x^^ ; (75a) 

gp{x) = exp - ^a;^) ; 9s' {x) = 6{x) . (75b) 
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Figure 1. A table which collects the most important quantities in discretized 
versions of QED, for regular lattices (right), and for arbitrary graphs (left). 



In our general frame of section^ this model corresponds to the matroid 

which is in the general form described in section |S1 Following the procedure described 
therein, we obtain for the dual matroid 



B 



I -J. 



S p 



pb 



in which sites are exchanged with dual sites, bonds are exchanged with plaquettes, and 
the relative orientation of bonds and plaquettes is reversed (in the geometrical duality 
for three-dimensional lattices, this corresponds to the fact that the orientations of 
bonds and dual bonds change in such a way that the vector product changes parity). 
Introducing the ad hoc notation for the partition function 

ZQED(e,i?)= /'d0,d^exp(-i^(Z?0)2-^^i;;2^ J]5((V.j%0, (76) 

•' b ^ p s' 

with the proper normalizations of the integrals, the duality relation states that 



■^QED(e,B) - / } ^QED [tT^B 

(47re2)— 



(77) 



where, as always, r is the rank of _B, and in our case is equal to the number of sites 
plus the number of bonds, and (n — r) is the rank of B, and in our case is equal to 
the number of dual sites plus the number of plaquettes. Again we see how duality 
exchanges a regime of high coupling with a regime of low coupling. 
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A concrete example could be the one of a cubic three-dimensional lattice, with 
the ordinary notions of sites, bonds and plaquettes, and the dual sites being the 
elementary cubic cells of the lattice. The geometric duality applied to the lattice 
naturally exchanges sites with dual sites and bonds with plaquettes. In this case, the 
presence of the delta function for the terms of current divergence is redundant, as the 
constraint is automatically satisfied. 

Variants in the spirit of disordered systems could involve a measure ^i{B) over an 
ensemble of random tessellations of the three-dimensional space, or, more generally, 
over an ensemble in which adjacency matrices are not directly derived from a 
geometrical structure, and the condition V • j = really needs to be enforced via 
a delta function in the Gibbs weight. 

In the case in which, following the recipe of section|Hl we have a self-dual measure 
fj,*{B), the fixed point of the involution for the electric charge e is at the value 2e^ = 1. 

14. Conclusions 

The notion of duality 'a la Kramers and Wannier, after its original introduction 
in 1941, has been applied to a wide range of periodic finite-dimensional systems of 
statistical mechanics. On the contrary, small work has been done in the direction of 
a systematic description of duality for a class of models sufficiently large to include, 
for example, mean-field disordered systems. 

Our idea in this paper is that, in the Kramers and Wannier formulation, a sort 
of two-step procedure is implicit: an algebraic duality, which exchanges variables and 
interactions, and performs a Fourier transform on the terms of the Gibbs weight, and 
a geometric duality, which, exploiting the geometrical structure of the lattice, allows 
for a reinterpretation of the dual system in a formulation which, in some contexts, is 
suitable for a direct comparison with the original system. 

In the more wide context of disordered systems the second step is in general 
impossible. Nonetheless, in many cases the first step is already suitable for a fruitful 
comparison of the two systems related by duality. 

The general recipe is the following: given some variables and some interactions, 
valued on some abelian group G, the Gibbs weight depends from two ingredients: 

• two sets of functions, weighting respectively the original variables, and the 
auxiliary variables corresponding to some linear combinations of the original ones 
(the interactions); 

• the specific pattern of linear combinations of the variables, encoded in a matrix 
J of group- homomorphisms. 

The dual system exchanges variable-functions and interaction-functions with the 
Fourier anti-transform and transform of interaction-functions and variable-functions, 
and the new pattern is given by the matrix —j"^ of transposed homomorphisms. 

We give some concrete applications for specific simple models: a random Gaussian 
model, a random Potts-like model, a random variant of discrete scalar QED. . . For 
each of these systems, a short discussion over the duality relation is proposed, and a 
hint over self-duality conditions is suggested. Actually, in the cases of interest, this 
discussion is a simple specialization of a general discussion done in section 8, and 
analogous to the one already present in ^ , where the duality, in its first application, 
was called variable-clause duality. 
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An unexpected connection with Matroid Theory has emerged in the analysis of 
the subject. The a postenori justification for this fact is that the pattern of interaction 
which is involved in the definition of the Gibbs weight is intrinsic with respect to a 
choice of basis for the space of configurations. Matroid Theory, in simple words, 
describes mathematical objects which convey basis-intrinsic information on matrix- 
like objects, like the pattern of linear dependence between the columns of a matrix. 
Hence, it is small surprise if the matroid related to the dual system is the matroid-dual 
of the matroid related to the original system. 

In a forthcoming paper we would discuss a further generalization of this 
framework, in which the crucial concepts of Fourier transform and of group- 
homomorphism for abelian groups are generalized to different integral transform and 
natural notions of homomorphism, applied to mathematical structures different from 
abelian groups. A first case is the one in which, for Partially Ordered Sets, the 
natural notion of transform is the Mobius transform, and the natural notion of 
homomorphisms is related to the theory of Galois connections |15) . 

Appendix A. Duality for graphical matroids and planar graphs 

An important subset of regular matroids| is the set of graphical matroids. Many results 
exist on the subject (cfr., for example, [H], chapt. 5), derived in the spirit of generality 
characteristic of Matroid Theory. In this appendix we will show some constructive 
procedures which apply to particular classes of graphical matroids, so in this sense 
they are less general, but more clear restatements of known facts. Hopefully, they will 
be useful to understand how to interpret the classical formulations of dualities for SM 
models on planar lattices inside our general setting of the duality. 

A precise terminology will be introduced in the body of the section. Here we 
resume the main results: 

Theorem 1 The "site-bond" matroid Jst and the "loop-bond" matroid Jib of a graph 
Q are dual matroids. 

Theorem 2 The "site-bond" matrix J' sb of a planar graph Q is equal to the "loop- 
bond" matrix J'lb of the corresponding geometrical dual graph Q* , and conversely the 
"loop-bond" matrix of Q is equal to the "site-bond" matrix of Q* . 

Appendix A.l. Matroid duality for graphical matroids 

Consider a graph Q, with vertices (or sites) s € V and edges (or bonds) b & E. We 
give an arbitrary orientation to the bonds. That is, for each bond 5, Sin(6) and Sout(^) 
will be respectively the first and second extremum of the bond. 



Define a 0-form ip as a function from V to a field F, and a 1-form -0 as a function from 
E to F. The operators d and d* transform 0-forms in 1-forms and vice versa with the 
following definition: 



s.„(b) 



b Sout (fc) 



(A.l) 



dip{b) = ip{sin{b)) - ifisoutib)) ; 



(A.2) 
(A.3) 




b 



X That is, matroids which are representablc over every field. 



General duality for abelian- group-valued statistical-mechanics models 



24 



Intuitively, they correspond to the discretization of the gradient operator on scalar 
fields and to the divergence operator on vector fields, generalized to arbitrary lattices. 

The sites arc implicitly a particular case of 0-forni, (/^(s) (s') = Ss_s' ■ The set of 
O-forms is a vector space over F, of dimension \V\. It is naturally decomposed in two 
vector spaces, <S and C, where O-forms in C are constant on each of the C connected 
components of G (so dim(C) = C), and O-forms in S arc obtained by quotient (so 
dim(iS) = S = \V\ — C). Projection on S is obtained via the operator ^d*d. 

Analogously, the set of 1-forms is a vector space B over F, of dimension B = \E\. 
The bonds are a particular case of generic 1-forms, where it is understood ip(^i,^(b') = 
Sbfii, and are a basis for the space B. We will define the loops as the 1-forms i/' such 
that d*^!) = 0. Elementary plaquettes are particular cases of loops: if dp is the border 
of an oriented plaquette p, ■)/'(p)(&) is the characteristic function on dp. The set of 
loops is a vector space £ over F, of dimension L. 

Note that S -\- L = B as & consequence of Euler formula. More deeply, the vector 
space B is decomposed into dS + C, that is, a generic 1-form t/j is uniquely decomposed 
into ip = d(p + X) with ip ^ S and x € £ (Hodge decomposition). 

The patterns of linear dependence of the set of bonds projected onto the vector 
spaces dS and C are matroids. They will be called respectively the site-bond and 
the loop-bond matroids. Given a choice of basis for S and £, we will denote the 
corresponding vectorial representations as Jgb and Jib- That is, for a set of O-forms 
{iPs} chosen as a basis for <S, the entries of Jst will be 

Jsb = dipsib), (A.4) 

and for a set of 1-forms {tpi} chosen as a basis for C, the entries of Jib will be 

Jib = Mb)- (A.5) 

We will give a particularly skill choice of basis, for which the two matroids are both 
in standard vectorial form, and manifestly dual. 

Choose a "ground" site ,so G V, and a spanning tree T C Q. Change the 
orientation of the bonds such that for each site s ^ sq there is exactly one out- 
going bond in E{T). Label this bond bg. Say that s s' iff an oriented path 7s'^s in 
T exists. Choose a basis for <S as 

Ms') = { (A.6) 

So we have 

{-1-1 S'^ Sinik), S i- Soutik) 
-1 sy Soutik), S i- Sinih) (A.7) 
o.w. 

We will choose a basis for C such that for each bond b E E \ E{T) there is exactly 
one 1-form ij^i such that 'ijji{b) is non-zero. So we can label the edges in \ E{T) as 
6;. More precisely, the 1-form corresponds to the only self- avoiding closed path on 
T U {6;}, oriented according to 6;. Analytically this can be expressed as 

S y Sin{bl), S i- Sout(&/) 
l\>l{bv) = \ l\)l{ps)={ +1 S ^ Sout(&z), S ^ Sin(60 (A.8) 

o.w. 
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It is immediate to check that Jsb and Jib are in standard vectorial form, respectively 
with the identity block on the left and on the right. From the observation that 
dips{bi) — —ipi(bs) we obtain that they are dual in matroid sense: 

Jsb = {I,Ysi); Jib = {-Yi^J). (A.9) 

As we have found an explicit choice of basis for which the standard form of the site- 
bond matroid and of the loop-bond matroid are dual, we have implicitly proved that 
the matroids are dual, which is the claim of Theorem 1. 

Appendix A. 2. Graphical duality for planar graphs 

In the previous section we have seen how the "abstract" site-bond and the loop-bond 
matroids for a given graph are dual in the sense of matroid duality. This is true in 
general, as the matroid definition is intrinsic, and specifically we have verified it in a 
particular choice of basis, for which both the matroids are in standard form. 

In this section we show the result of Theorem 2, that the site-bond and the loop- 
bond incidence matrices for a connected planar graph and for its graphical dual are 
exchanged by the graphical duality. This result relies on a particular choice for the 
loops as the elementary plaquettes of the lattice, which is the more natural for the 
case of planar graphs, but in this sense is more accidental. 

A graph Q is said to be planar when it is embeddable on a genus-0 surface with 
no crossings between the bonds. Consider the case in which Q is connected. The 
classical Euler relation for this kind of graphs claims that V + F = B + 2, where V is 
the number of sites, B is the number of bonds and F is the number of faces. 

One can immediately define the site-bond and loop-bond incidence matrices, Agb 
and Alb, for non-oriented bonds and loops. Just state that Agb is 1 if the site s is an 
extremum of the bond b and elsewhere, and Aib is 1 if the bond 6 is a side of the 
face I and elsewhere. 

The graphical dual Q* is defined by the following procedure. Replace each 
plaquette I by a dual vertex s* (white bullets in Fig.^, each bond 6 by a dual (dashed) 
bond b*, which crosses it and connects the dual vertices of the two neighbouring 
plaquettes, and each site s with the dual plaquette 1* delimited by the dual bonds of 
the bonds outgoing from the site. It is clear that 

As'b'{g*) = Aib{g); Ai,b'{g*)^A,b{g)- (a.io) 

If we put an orientation on the bonds and a direction of rotation on the faces, the 
matrices Asb and Aib will be promoted to some matrices J^^ and J/^, in which some 
minus signs appear. Again because of the planarity of the graph, we could give a 
recipe such that the relation (|A.10p above is promoted to a relation 

j's'b'igi^j'ibig); J'l'b'ig*) = j'sbiQ) ■ (a.h) 

The recipe is the following: choose the direction of rotation to be the same for all the 
faces, (for example clockwise), and opposite between faces and dual faces. Choose the 
orientation of the bonds arbitrarily. The orientation of the dual bonds must be such 
that the vector product between each bond and its dual bond has the same sign with 
respect to the normal versor outgoing from the surface. 

Both the matrices J^^ and J^'^ are not of maximal rank, and the only linear relation 
is given by the sum of all the rows with equal coefficients. This accounts for the offset 
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Figure Al. Graphical dual for a planar graph (left), and a consistent choice of 
orientations for the bonds and directions of rotations for the faces for the relation 
lA.lH to hold (right). 

of 2 in the Euler formula V F = B -\- 2, in fact the number of "independent sites" 
in the loop-bond matroid (0- forms up to an arbitrary constant) is S" = — 1, and 
the number of independent loops is L ~ F — 1, from which we recover the relation 

S + L = B. 
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